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Abstract 



We discuss a generalized quantum microcanonical ensemble. It 
describes isolated systems that are not necessarily in an eigenstate of 
the Hamilton operator. Statistical averages are obtained by a combi- 
nation of a time average and a maximum entropy argument to resolve 
the lack of knowledge about initial conditions. As a result, statisti- 
cal averages of linear observables coincide with values obtained in the 
canonical ensemble. Non-canonical averages can be obtained by tak- 
ing into account conserved quantities which are non-linear functions 
of the microstate. 

1 Introduction 

In a recent paper, Goldstein et al [T] argue that the equilibrium state of a 
finite quantum system should not be described by a density operator p but 
by a probability distribution f(ifj) of wavefunctions ip in Hilbert space 7i. 
The relation between both is then 




(1) 



Integrations over normalised wavefunctions may seem weird, but have been 
used before in [2] and in papers quoted there. Goldstein et al develop their 



arguments for a canonical ensemble at inverse temperature (3. For the micro- 
canonical ensemble at energy E they follow the old works of Schrodinger and 
Bloch and define the equilibrium distribution as the uniform distribution on 
a subset of wavefunctions which are linear combinations of eigenf unctions of 
the Hamiltonian, with corresponding eigenvalues in the range [E — e, E + e]. 
The parameter e is fixed and small. However, there is no compelling reason 
why an isolated system should be in an eigenstate of the Hamiltonian, or in 
a superposition of stationary states all with about the same energy. On the 
contrary, time-dependent states are generic — their time average determines 
the equilibrium state. Such a time-dependent state cannot be decomposed 
generically into eigenstates with energy in the range [E — e, E + e] . Its wave- 
function only satisfies the conditions 

ifj) = 1 and Hifj) = E. (2) 

Hence, in a generalized microcanonical ensemble, it is obvious to define the 
microcanonical density operator p by means of where the integration 
extends over all normalised wavefunctions ip satisfying (j2J). This approach 
can be worked out explicitly, as has been done in a recent paper by Brody 
et al [3 E] . However, we believe that this is not yet the correct way of defin- 
ing the generalized quantum microcanonical ensemble. Indeed, the density 
operator p obtained in this way has non-physical properties, even for simple 
rz-level systems. E.g., the heat capacity of the 2-level system does not depend 
on energy, contradicting experimental evidence that heat capacity vanishes 
as the energy approaches the ground state energy. The 4-level system ex- 
hibits a phase transition, which is totally unexpected. In our opinion, this 
approach has to be modified in two aspects: 1) Due to the occurrence of 
non-commuting observables replacing time averages by phase space averages 
is not straightforward; 2) The abundance of conserved quantities in quantum 
systems must be taken into account. 

Recently, the interest in the microcanonical ensemble has increased be- 
cause of the possibility of microcanonical phase transitions, occurring even 
in finite systems jH]. The phenomenon of negative specific heat was known 
in astrophysics. Its existence has been demonstrated rigorously in the model 
of Hertel and Thirring [B]. When such an instability occurs then the usual 
argument of equivalence of ensembles becomes invalid. Hence, one cannot 
anymore profit from the freedom of choice to use the most simple ensem- 
ble, but one is forced to do microcanonical calculations. Up to now these 
are almost exclusively done using classical mechanics. However, there are 
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problems, such as the stability of metal clusters |H] or of atomic nuclei [H] , 
that one would like to treat quantum mechanically. The standard definition 
of the quantum microcanonical ensemble is useless for this purpose because 
it limits itself to (nearly) stationary states. By including time- dependent 
states the behaviour of the system approaches that of its classical equivalent. 
This raises the hope that a microcanonical phase transition, present in the 
classical system, shows up in the quantum system as well. 

Conserved quantities are discussed in the next section. In Section 3 mi- 
crocanonical averages are defined. The classical ergodic theorem is used to 
replace time averages by integrations over wavefunctions. In Section 4, the 
non-uniqueness of the equilibrium state is resolved by means of the maxi- 
mum entropy principle. As a result, the microcanonical averages coincide 
with canonical averages. This takes away the need of proving the physical 
relevance of the present approach. Indeed, the evidence that the canonical 
ensemble gives an adequate description of physical phenomena is overwhelm- 
ing. One should therefore expect only minor differences with the standard 
treatment of quantum statistical mechanics. One of these is that the micro- 
canonical ensemble is described by a set of equiprobable wavefunctions, and 
not only by a density operator. A new feature is the possibility of calculating 
statistical averages of nonlinear functions of the state of the system, while 
quantum expectations of a quantum observable always depend linearly on 
the state of the system. 

Sections 5 and 6 analyse the role of additional conserved quantities for 
the example of the harmonic oscillator. Coherent states, modified with loga- 
rithmic corrections, are introduced. The final section contains a preliminary 
evaluation of the present approach. 

2 Time averages and lack of ergodicity 

The ergodic theorem plays an important role in understanding the micro- 
canonical ensemble for classical (i.e. non-quantum) systems. Even when the 
ergodic hypothesis is not satisfied the statement remains that first averag- 
ing over time and then over initial conditions in phase space gives the same 
result as directly averaging over phase space alone. The only condition for 
this to be true is that the probability distribution over phase space is time 
invariant. In the quantum context the time average cannot be omitted and 
is an essential step in obtaining a density matrix which commutes with the 
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Hamiltonian. 

Assume a Hamiltonian H with energy eigenvalues E n and corresponding 
eigenstates ip n . Let ip = J2 n ^n^n be a wavefunction of the desired energy 

^W) = Y,ME n = E. (3) 

n 

Its time evolution is given by 

It = Ke~ ih ~ lEnt ^ n . (4) 

n 

Then the time averaged expectation value of an operator A is given by 

lim i F dt (MArk) = Tr p+A (5) 

with 

Pip = Y\ :K n\ 2 \4>n)(4>n\- (6) 
n 

With a strict definition of (generalised) microcanonical ensemble, any density 
operator of the form (0), with A n satisfying ^nl^™| 2 = 1 an d ©, is a 
candidate for describing the equilibrium distribution. But this is clearly any 
density operator p which is diagonal together with H and which satisfies 
Tr pH = E. This abundance of stationary states is a consequence of the 
large number of conserved quantities. From a classical point of view this 
means that a quantum system is almost always non-ergodic. 

3 Microcanonical averages 

An example of a non-linear function of the state of the system is 

f{^3) = mABip) - (VW> <</W>, (7) 

where A and B are quantum observables. The statistical average of such a 
function is defined by 

(/>*= lim \ Fdtffatft), (8) 
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whenever this limit converges. Note that if) = J2 n An^n with Hif> n = E n ip n 
implies that if> t = J2 n exp(— ifir 1 E n t)if> n . Hence, by the classical ergodic 
theorem one has 



Jd X f (^Xne^^Xne-^if)^ , 



(fh= I d X .f[ ) j A ! />v, i .)^V^ rj; j . («)) 
where 



II 



The latter is an integration over wavefunctions, as in However, the 

domain of integration is limited to a time-invariant subset of wavefunctions 
generated by a single if). Note that the ergodic hypothesis, needed for the 
validity of 0, is satisfied in the generic case. If it is not satisfied then, in line 
with the tradition of classical statistical physics, Q is taken as the definition 
of the microcanonical average rather than (JSJ). 

In the linear case is f(if>,if>) = (if>\Aif>). Then (jHJ) becomes 



j d X ^X m e^i) m \A^X n e^i) n \ 

\ m n I 

= ^\X n \ 2 (^ n \Aip n ) 

n 

= Trp^A. (11) 
The average of the non-linear function (JZj) becomes 

Tr f^AB- I d X (^X m e i ^if) m \Aj2^ne iXn 'ipn) 
X m ^ Xm 'ifj m ^B^ X n <e ix "'ip n , \ 

\ m 1 n' I 

= Tip^AB - (Tr pi,A)(Tr pi,B) 

-Tr H A H B + \X n \ 2 (^n\Aip n ) (ipn\Bif) n ). (12) 

n 

The last two terms are a consequence of the non-linearity. 
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4 Maximum entropy 



The non-uniqueness of the micro canonical equilibrium state can be lifted in 
more than one way. One approach would be to integrate (0) over all nor- 
malised wavefunctions ip satisfying @. This can be worked out but leads to 
some strange results. In particular, the resulting density operator p depends 
non-analytically on E. We do not pursue this approach. 

The lack of information about the coefficients A n can be dealt with by 
means of the maximum entropy principle. Let p n = |A n | 2 . Take entropy of 
the form 



as usual. After introduction of the Lagrange multipliers a and (3, to take the 
constraints ^2 n p n = 1 and ^2 n fJ> n E n = E into account, one obtains 



But then the density operator p^, defined by (JHJ), coincides with the Boltzmann- 
Gibbs-von Neumann result of the canonical ensemble. 

This choice of microcanonical wavefunctions has many advantages. First 
of all, most microcanonical results coincide with the canonical ones. This 
eliminates the need to convince the community that the microcanonical re- 
sults are physically acceptable. Further advantages follow from the use of 
the maximum entropy principle. The canonical density operator is the most 
likely density operator satisfying the energy requirement. It minimises the 
amount of information produced by the model. But in addition, it guarantees 
thermodynamical stability. Indeed, the entropy (fT3~j) can be identified with 
thermodynamic entropy. Its derivative with respect to energy E is inverse 
temperature, and equals k B [3. 

The choice ()14|) appears to be very satisfactory. Nevertheless, some fur- 
ther refinement may be required, as will be made clear in the next sections 
for the example of the harmonic oscillator. 

5 The harmonic oscillator 

We now start to develop an alternative approach based on the example of 
the harmonic oscillator. The next section shows then how to reconcile this 




(13) 



n 



(14) 
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approach with the maximum entropy approach discussed above. The Hamil- 
tonian of the harmonic oscillator is 

H = ^-P 2 + \mu?Q 2 (15) 
2m 2 

Its eigenvalues are E n = hu>(^ + n), n — 0, 1, • • •. 

In the spirit of Boltzmann's entropy we look for a measure of the phase 
space extent of the set of wavefunctions ip t ,t G R. To this purpose, consider 
a semi-classical phase portret obtained by measuring average position Q and 
average momentum F. Note that the exact solution of Heisenberg's equations 
of motion is known 

Q t = Q cos ujt-\ — — Fsinu;t (16) 
mtu 

P t = —muQ sirnot + P cos cut. (17) 

Hence, for any if), the averages (Q) t = (ipt\Qi>t) — (^iQt^) and (P) t = 
{ip t \Pipt) = {ip\Ptip) describe ellipses in the ((F), (Q))-plane. The logarithm 
of the 'size' of the ellipse can be taken as definition of (non-equilibrium) 
entropy. It is well-known that the ellipse is maximal when ip is a coherent 
state wavefunction, i.e. 

^ = e-l^l 2 / 2 -r=z^ n (18) 
n=o vn! 

for some complex z. The maximal values, that are attained, equal 



2h 

(<5)max = Aq|z| With Xq = \l (19) 

V muo 

(P)max = Ap|z| with Ap = \ / 2hmuj. (20) 
The energy condition for such a wavefunction reads 

F = ^Q + M 2 ). (21) 

It is now obvious to claim that the coherent states are the equilibrium 
states of the microcanonical ensemble. But this claim is incompatible with 
the maximum entropy argument of Section 4. Indeed, it implies a Poisson 
distribution of the amplitudes |A n | 2 instead of an exponential distribution. 
The way out of this conflict is the introduction of a microcanonical ensemble 
with two conserved quantities. 
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6 Additional conserved quantities 

Note that classical energy 



(22) 



is a conserved quantity of the quantum harmonic oscillator. Hence it is mean- 
ingful to study a microcanonical ensemble where both E c \(ip) and quantum 
energy E = (^\H^) are fixed. The equilibrium state of this ensemble is 
described by the wavefunction if; = X n ifj n obtained by maximising 



- |A n | 2 ln |A n | 2 - aJ2 l^| 2 - PJ2 \X n \ 2 E n - yE^). 



(23) 



Introduce creation and annihilation operators a and a*, as usual. A short 
calculation shows 



with 



£7 d (V>) = hu(il>\a*il>) (if>\ail)) 
= fuv\C\\ 



C = ^2 ^n+iKVn + 1. 



(24) 



(25) 



Variation of (j23j) with respect to |A n | then gives 
= -ln\\ n \ 2 -l^a- [3E n 

V Xn X r 



(26) 



This equation can be read as a recurrence relation for the coefficients A n . 
Fixing Ao, it determines a unique normalised wavefunction, which is denoted 
ijj(f3,j, Aq). Let us analyse the asymptotic behaviour of the A n . Put 



A, 



-j(huj — 



(27) 



with x n real and with c n = Ylp=o(VP^ n y/P)- Then fl2^j) becomes 
./■„ 2 In c n — 2 In x n — rain ( -^{(^(hu) 2 J — 1 — a — (3hw 



(H 
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(28) 



From this expression follows that the x n are constant for large n (note that 
21nc n increases as rzlnn). Hence, (12 7|) shows that the A n behave asymptot- 
ically in the same manner as the expansion coefficients of a coherent wave- 
function, be it with logarithmic corrections because c n 7^ \fn).. 

7 Discussion 

We started with the observation that one cannot simply replace quantum 
time averages by averages over some phase space of quantum mechanical 
microstates. Indeed, time averages are essential to obtain a statistical de- 
scription by a density operator which commutes with the Hamiltonian of 
the system. But, restricting statistical mechanics to a formalism based on 
nothing but time averages is not very satisfactory because of the intrinsic 
non-ergodicity of quantum systems. The non-unicity may be resolved by 
means of the maximum entropy principle. As a result, the microcanonical 
ensemble is described by the same density operator as the canonical ensem- 
ble. This is not bad in itself because the results of the quantum canonical 
ensemble are widely believed to agree with experimental reality. However, 
this cannot be the whole story. We expect that microcanonical systems can 
be thermodynamically unstable, while canonically this can only happen in 
the thermodynamic limit. 

The microcanonical ensemble does not any longer coincide with the canon- 
ical ensemble when constraints are introduced, which depend non-linearly on 
the microstate. We have tried out this idea on the harmonic oscillator. Based 
on semi-classical arguments one expects that the coherent states are the mi- 
crocanonical equilibrium states of the harmonic oscillator. Hence the obvious 
question is wether there exist conserved quantities such that coherent states 
are obtained by maximising entropy under the constraint that the conserved 
quantities have given values. The quantum harmonic oscillator has a con- 
served quantity which, for convenience, we have called the classical energy. 
It controls whether the microstate is more quantum-like or more classical. It 
vanishes for eigenstates and has its maximal value for coherent states. For 
intermediate values, the microcanonical equilibrium states behave asymptot- 
ically in the same way as coherent states. 
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The micro canonical ensemble describes closed systems. This means that 
the system plus environment is in a product state. This state must however 
be prepared by the experimenter starting from a situation where the system 
interacts with its environment. It has been shown recently that under rather 
general conditions on system plus environment the reduced density operator 
of the system is that of the canonical ensemble [1U| . It is then no surprise 
that, after isolating the system, it is still described by a density operator of 
the canonical ensemble. The only correction to this view that we want to 
make is that one should take into account additional conserved quantities, 
for example, those which control the degree of quantumness of the system. 
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